function dt_new = varTimeStep(V_new, V_old, dt_old, D, dnorm)
%VARTIMESTEP Summary of this function goes here
%   Calculates new dTao using variable timestep selector

% Find MaxRelChange
RelChange = abs(V_new - V_old)./max(max(D, abs(V_new)), abs(V_old));
MaxRelChange = max(RelChange);
% New dTao
dt_new = dnorm/MaxRelChange*dt_old;

end

